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Abstract 

The proposal that the interaction between a macroscopic body and 
its environment plays a crucial role in producing the correct classical 
limit in the Bohm interpretation of quantum mechanics is investi- 
gated, in the context of a model of quantum Brownian motion. It is 
well known that one of the effects of the interaction is to produce an 
extremely rapid approximate diagonalisation of the reduced density 
matrix in the position representation. This effect is, by itself, insuffi- 
cient to produce generically quasi-classical behaviour of the Bohmian 
trajectory. However, it is shown that, if the system particle is initially 
in an approximate energy eigenstate, then there is a tendency for the 
Bohmian trajectory to become approximately classical on a rather 
longer time-scale. The relationship between this phenomenon and 
the behaviour of the Wigncr function post-decoherence (as analysed 
by Halliwell and Zoupas) is discussed. It is also suggested that the 
phenomenon may be related to the storage of information about the 
trajectory in the environment, and that it may therefore be a general 
feature of every situation in which such environmental monitoring 
occurs. 
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1. Introduction 

It is well-known |Q, ||, ||, ^ that the trajectories in the Bohm Interpretation of 
Quantum Mechanics are often highly non-classical. This gives rise to an important 
problem for the Bohm interpretation: namely, the question as to how the inter- 
pretation can account for existence of generically quasi-classical trajectories on the 
macroscopic level. In a previous paper Q we argued that the Bohm interpretation 
typically fails to produce the correct classical limit if the system is isolated. The 
purpose of this paper is to argue that the correct classical limit does emerge once 
one takes into the account the effect of the environment. Our discussion extends 
the analysis given in Chapter 8 of the book by Bohm and Hiley Q . 

We are particularly interested in the connection between the role of the envi- 
ronment in the Bohm interpretation, and the phenomenon of decoherence, which 
plays a central role in the decoherent histories approach [f[ |[ 0, || , and in Zurek's 
existential interpretation H, The mechanism considered by Bohm and 

Hiley — the scattering of a beam of radiation or other particles — is also one of the 
mechanisms by which decoherence is produced |§, ^ This has suggested to 
some authors |l^ that the process by which the Bohmian trajectory becomes 
quasi-classical is closely related to the phenomenon of decoherence. The suggestion 
is certainly plausible. However, it is not entirely clear, just from the argument 
given by Bohm and Hiley, that the suggestion is actually correct. At the macro- 
scopic level decoherence is an ubiquitous phenomenon, which can be produced by 
a wide variety of different mechanisms. By contrast, Bohm and Hiley only con- 
sider the particular case of a scattering process. One would like to know whether 
other kinds of interaction between a macroscopic body and its environment also 
have the effect of causing the Bohmian trajectory to become approximately classi- 
cal. More generally, one might ask whether every process which causes decoherence 
also causes the Bohmian trajectory to become approximately classical, or whether 
it is only some of them. These are the questions which motivated our investigation. 

We will focus on the models of quantum Brownian motion which have been 
discussed by Caldeira and Leggett, Hu, Paz and Zhang, and many others 
p^ . These models have played an important role in studies of decoherence, and are 
therefore a natural starting point for an investigation into the role of decoherence in 
the Bohm interpretation. We accordingly consider a system particle, with position 
X and momentum p, interacting with a heat bath consisting of N other particles 
with positions Xi, . . . , xn and momenta pi, . . . ,pn- The Hamiltonian is 

^ = (i^' + I'^'^oi') + E + ^mrivlxA + E ^rXXr (1) 

\ / r—1 ^ ^ r— 1 

Here uq denotes the bare frequency of the system particle. The renormalised fre- 
quency will be denoted uj. 

The model is characterized by the spectral density 
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I{io')^^-^S{u;'~u;r) (2) 

r—1 

Taking / oc w' (for uj' < the cut-off frequency) gives the Caldeira-Leggett model [|6| . 
If one leaves / arbitrary one obtains the general class of master equations derived 
by Hu, Paz and Zhang 

At i = the heat bath is taken to be in the thermal state with density matrix 

Pbath = J^cxp I j (3) 
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where 

N 



r—1 ^ ^ 



and is a constant. We assume that at i = system+heat bath are in the product 
state |'0sys) (V'sysl i^ipbath- At i > the density matrix describing system+heat bath 
wiU consequently be 

In the conventional approach one now integrates out the environmental degrees of 
freedom, and focusses on the behaviour of the reduced density matrix. Unfortu- 
nately matters are not so simple in the Bohm interpretation. 

In the Bohm interpretation a mixed state such as phath is taken to describe an 
ensemble of pure states 

/Sbath = If/lQ) ((/IqI (4) 

a 

This way of writing the density matrix is not simply a mathematical device, as in 
the conventional approach. Rather, one takes it that at i = the heat bath actually 
is in one of the pure states in the ensemble, with pa being the probability that it 
is in the state I^^q). The problem we then face is, that the density matrix does not 
uniquely determine the ensemble, and that in the Bohm interpretation it makes a 
difference which ensemble we choose (for a classification of the set of all discrete 
ensembles corresponding to a given density matrix see Hughston et al ||2^). We 
discuss this point further in Section ^. 

Suppose that a particular ensemble has been chosen, and suppose that at t = 
the heat bath is in the pure state \(t>a)- Then at t > system+heat bath will be in 
the pure state 

=6-"*^/" l^sys)® (5) 

and the Bohmian velocity of the system particle will be given by 

(a),, ^ _ ^Im((^'a(t) I X,Xi, ...,Xn)-§^{x,Xi,...,Xn\ 

(r, x, xi, . . . , xatj — ■ "2 

m|(a;,a;i, ...,xn \ ^a{t))\ 

(6) 

We see that the Bohmian velocity depends, not only on x, but also on . . . , a; at, 
as well as the index a. The reduced density matrix clearly does not provide enough 
information to calculate this function. Consequently, the problem of determining 
the effect of the interaction with the heat bath is significantly more difficult in the 
Bohm interpretation than it is in the conventional approach. 

Nevertheless, although the reduced density matrix does not provide us with 
complete information regarding the Bohmian velocity of the system particle, it 
does tell us something. To see this, consider the effect of averaging over all the 
possible values of xi, . . . , xn, and of the index a: 



VEit,x) = 



J2a Pajdxi... dXN 


\{x,xi, . . . 


, Xn 




2 

V 


B 


Xi,. . . ,Xn) 


J2a Pa I dxi . . . dXN 


\{x,xi, . . . ,xn\ 
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We will refer to ve as the ensemble-averaged velocity. The reduced density matrix 
elements are given by 

(a;| prcd{t) \x') 

= ^Pq / dxi . . .dxN {x,Xi, . . . ,XN \x' , Xi, . . . , Xn) 

It is then straightforward to infer 

- t, . ft Im(^(x| p^t) \^')l,=J 
= m {x\ p^^t) \x) 

from which we see that the reduced density matrix does provide enough information 
to calculate we- 

In this paper we will investigate the behaviour of we and Wg"'' as functions of 
time for the case when the initial system state is an approximate energy eigenstate, 
of the form Q 

li^sys) = ^ Cr\n + r) (8) 

In this expression |n) denotes the n*'' eigenstate of the isolated system particle 
Hamiltonian: 

Hsys \n) = En \n) 

where iJgys — P^/ (2to) + moj'^x^ /2 and £"„ = [n + 1/2) huj. We assume that n':^ 1 
(so that the state is highly excited), and An ^ n (so that the energy distribution 
is sharply peaked about the mean). Classically one would therefore expect the 
particle to be following a well-defined orbit with energy close to and amplitude 

1/2 

close to Xmax = {2En / {mijo'^)) . In particular, when the particle is located at x, 
one would classically expect its velocity to be close to ±pci{x)/m, where 

Pciix) ^muj (x^ax-a;^)' 

On the other hand it was shown in ref. [Q that, assuming the system to be isolated, 
there is only a high probability of this being true of the Bohmian velocity at all 
stages of the motion in the very special case for which |'0sys) is a narrowly localized 
wave packet. In the following we will show that the effect of the interaction with 
the heat bath is to make the distribution of Bohmian velocities eventually become 
approximately classical, whether or not this is true initially. 

We begin, in Sections ^ and ^, by considering the behaviour of the ensemble- 
averaged velocity we- The feature of the interaction with the environment which has 
probably attracted most attention is the tendency of the reduced density matrix to 
become approximately diagonal in the position representation. In Section ^ we show 
that this phenomenon is, by itself, insufhcient to produce approximately classical 
behaviour of the Bohmian trajectory. However, the interaction has other important 
effects, apart from the approximate diagonalisation of prcd- In particular, Halliwell 
and Zoupas ||2l| have shown that, in the case of the Caldeira-Leggett model, the 
Wigner function becomes non-negative after a sufficient elapse of time- In Section ^ 
we show that, as a consequence of this effect, we comes to lie approximately within 
the classical range 

-Pciix)/m < veIx) < pc\{x)/m (9) 

We also derive conditions for this to occur in the case of other models of the type 
defined by Eq. (|l|) (and, in fact, for a number of models which are not of this type). 
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Inequalities (^) represent a necessary condition for the Bohmian trajectory to 
be approximately classical. However, they are clearly not sufficient. In Section ^ 
we accordingly calculate the function Wg"'' (t, ) on the assumption that 

the ensemble described by pbath consists of coherent states [see remarks following 
Eq. (U)]. Our calculation is based on Halliwell and Yu's alternative derivation 
of the Hu-Paz-Zhang master equation (also see Anglin and Habib |^), which has 
the advantage (from our point of view) that, unlike the usual path integral meth- 
ods, it allows us explicitly to keep track of the heat bath degrees of freedom. We 
show that in the case of the Caldeira-Leggett model, for sufficiently large values 
of t, there is a high probability that Vq will be close to one of the classical values 
±pci{t, X, xi, . . . , xj^). We also derive conditions for this to occur in the case of 
other models of the type defined by Eq. (|l|) 

Finally, in the conclusion, we discuss the bearing that these results have on 
the questions which provided the original motivation for this investigation, and we 
suggest some directions for further enquiry. 

2. Effect of Approximately Diagonalising p^cd in the 

X- REPRESENTATION. 

One of the most striking effects of the interaction between a macroscopic body 
and its environment is that the reduced density matrix tends rapidly to become 
approximately diagonal in the position representation. We begin by showing that 
this effect is not, by itself, sufficient to cause the Bohmian trajectory to become 
quasi-classical. 

The point is most conveniently illustrated in the context of the Caldeira-Leggett 
model, for which the master equation is [ p6| 

ih—prcd = [Hsys,Prcd\ +7[X, |p,prcd|j [x, [x , Prcd\\ (10) 

where i?sys is the renormalised system particle Hamiltonian (which, for the purposes 
of this section, need not be assumed to be of oscillator form), and where {.,.} 
denotes an anti-commutator. It should be noted that this equation is not exact. 



and that it does in fact violate the positivity of picd over very short times 24, p5|. 
However, it provides a good approximation over somewhat longer times. 

Under conditions where the last term on the right hand side of Eq. ( p^ ) dom- 
inates, and provided that t sufficiently small [but not so small as to render the 
approximation of Eq. ( p^ ) invalid] , one has [^ |l^ 

Predit) W) « exp [~Atix ~ X'y] {x\ p,ed(0) \x') (11) 



where A = {2m'^k'QT) j is the localization rate. In the case of a macroscopic object 
A is typically very large [|l2|, p^, even when the interaction with the environment 
is comparatively weak. Eq. ( |l l| ) consequently plays an important role in attempts 
to explain the emergence of an effectively classical statistics of "facts" [ 27| f rom an 
underlying theory which is fully quantum mechanical (|, 0, | |[ ^ p^, |13|. 



Substituting the expression given by Eq. (|11|) into Eq. ([7|) wc find 

VE(i,a;) w «E(0,a;) (12) 

However, at t = system-|-environment are in the product state [V'sys) (V'sysl ®/5bath, 
which means that we(0, x) coincides with the actual Bohmian velocity of the system 
particle a,t t — 0, 

h Im((V'sys \x)^{x\ V^sys)) 

M^.x) = ^ — ■ ^ ^ 
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It was shown in ref. that, for many choices of |'0sys), this quantity tends to take 
values greatly in excess of the classical speed. In view of Eq. ( [l^ ) the same must 
be true of VEit,x). It follows, that the approximate diagonalisation of prcd in the 
x-representation is not, by itself, sufficient to produce generically quasi-classical 
behaviour of the Bohmian trajectory. 



3. The Behaviour of ve at Later Times 



The approximation of Eq. (11) is only valid for sufficiently small values oft. We 
now want to investigate the behaviour of the Bohmian velocity over longer time- 
scales, and for other models of Brownian motion, apart from the Caldeira-Leggett 
model. We will consider the function VE{t,x) in this Section, and the function 

iig""* (t, X, xi, . . . , Xn) in Section ^ 

The advantage of considering v-Eit,x) is that in order to calculate it one only 
needs to know the master equation. We will consider master equations of the form 

d 1 

iS^Prod = -^[{hl{t)x'^ + h2{t)p^ + hy,{t) {x,p}) , /5,.cd] + l{t)[x, {p,/5red}] 

i 

" ft' 

for which the right-hand side is quadratic in x and -p. This class includes equations 
of the Hu-Paz-Zhang jl^ type, corresponding to Brownian motion models of the 
kind defined by Eq.(|^). It also includes those equations of the Lindblad form 
for which the right-hand side is quadratic in x and p. In particular, it includes the 
equation discussed by Diosi ||2^, |2|] . 

It is most convenient to work in terms of the reduced Wigner function. 



■{j\\{t)\f), [p, prcd]] -2Ji2(t)[i, [p,pi.cd]] + 'hi{i)\i, [x,/5,.cd]]) (13) 



Wrod(i,a;,p) = \ I dy exp[ ^py 



y 

X 

2 



/5rod(t) 



Expressing Eq. (y_3|) in terms of IVrcd we find 

2 



dt 
where 



(14) 



(15) 



»7 = 
K(t) = 

m = 



'Jli(t) Ji2(t)^ 
yi2(t) J22(t), 

It is straightforward to verify that the solution to Eq. ([l5[) may be written M, 09, |3C 



1 



Trdet A(t)^detM(t) 



dVexp ~ [t]' - A-\t)'nf M-\t) {q' - A-\t)'n) l^rcd(0,r,') (16) 



where the matrices A and M are defined by the equations 

|-M = 4A-ij(A-i)^ 
dt ^ ' 



(17) 
(18) 
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(superscript "T" signifying "transpose") subject to the initial conditions 

A(0) = 1 M(0) = 

Before proceeding further it wiU be useful to relate this equation to the discussion 
in the last section. Specialising to the case of the Caldeira-Leggett master equation, 
Eq. (p^), with H taking the oscillator form j? I(2jn) + mw^a;^/2, one has 

/ -i^ 

[ m 


D 

where D = 2m-ik^T . If wt, 7i < 1 Eqs. and (|l|) then imply 

1 

^mcj^t 1 



1 

2m 



(19) 



so that 



Wrod(i, a;,_p) 



dx'dp' exp 
3m / 



Dt^ ^ 



Wrcd(0,a;',p') (20) 



where we have made the further approximation A^^(<) w 1. In this expression the 
width of the Gaussian convolution in the p-direction is cx , whereas the width 
in the x-direction is oc ^2. It follows that, if t is sufficiently small, there will be a 
significant degree of smoothing in the p-direction, but no significant smoothing in 
the x-direction. To be specific, suppose that the initial system state IV'sys) is of the 
form specified by Eq. (ph, and suppose that t <$^tc, where 



tr. 



D 



(21) 



where Ab — h/ {mLox^sx) is the minimum value of the de Broglie wavelength. In 
that case Wjod will be nearly constant over the width of the Gaussian in the x- 
direction, and we can approximately write 



W^f,d{t,x,p) 



J dp' i J dx' exp 



3m^ 



(x'-xY 



3m 



-J^i^'-"") {P' -P) 
( 1 



cxp 



^11 \2 



I^rcd(0,X,p') 



V IvrZJt 



dp' exp 



{^-pf 



W,ed(0,X,p') 



If this result is re-expressed in terms of the density matrix one recovers Eq. ( |11| ) 
(with h. = D/ll?). It follows that the discussion in the last section only applies to the 
situation when i <C tc, before there has been any significant degree of smoothing 
in the cc-direction. The question we now have to consider is whether there is a 
tendency for the distribution of Bohmian velocities to become generically quasi- 
classical when t > tc- 

The result established by Halliwell and Zoupas |^ provides some preliminary 
indication that such an outcome might be expected. Halliwell and Zoupas show 
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that, in the case of the Caldeira-Leggett model as apphed to a free particle, with 
negligible dissipation, the Wigner function becomes strictly non-negative once t > 
(3/16)^^^ iioc, where t\oc is the localization time given by 

1 



(22) 



It is easily seen that this is also true for the case of a harmonically bound particle 
considered here (provided cj^ioc, T^ioc <C 1, so that the approximation of Eq. ( pO| ) 
is still valid at t = tioc)- More generally, what is essentially the same argument 
shows that, whenever the Wigner function propagator takes the form specified by 
Eq. dl^), Wrod becomes strictly non- negative once 

detM(i) > 

It follows from Eqs. (0) and (|l^ that 

/ dppWrcdit,x,p) 



VEit,x) 



(23) 



TO / dpWrcd{t,X,p) 

If Wred is non-negative, and if it is negligible outside the region enclosed by the 
classical energy surface p = ±pci{x), then it can be seen from this expression that 
ve will lie within, or close to the classically permitted range for the mean velocity: 
-Pci{x)/m < VEit, x) < pc\{x)/m. 

We now illustrate this phenomenon by calculating Wiod when IV'sys) is an ap- 
proximate energy eigenstate of the form specified by Eq. (||). It turns out that for 
such states we typically comes to lie within the classically permitted range of values 
at times significantly earlier than t\oc- 

It was shown in ref. B, using the WKB approximation, that 



{x I V'sys) ~ i exp 



-Six) 



9-{x) - exp 



:S{x) 



provided that x is not close to one of the classical turning points at a; = ±Xmax- In 
this expression 

h 



Six) 



dx' pc\{x') -f 



and 

9±{x) 



(27rp,i(a;)) 



A„ Cr exp 



±ir sin 







if a^max X <i a^max 

otherwise 



This gives, for the reduced Wigner function at i = 0, 
VFrod(0,a;,p) 



\jdy e^pUpy 



X •j^exp 
X < exp 



exp 



-S 



Substituting this expression into Eq. ([l^) and carrying out the p'-integration gives 
W,cd{t,x,p) 



X exp 



dx'dy 



X ■j^exp 
X I cxp 



.9- 1 2^' " |- ) - exp 



-5 X + — 



g_ [ X + y ) - exp 



, y' 



-5 X + — 



9+\x - — 



* I / , y 



(24) 



where we have assumed that t is sufficiently smah to justify the approximation 
A^^{t) sa 1, and where we have set 



detM-i = A 



(25) 



In order to evaluate this expression we note, first of all, that the functions g± are 
slowly- varying Q . We may therefore write 



9-{x' ~ 9-{x) 



provided that the Gaussian peaks are sufficiently narrow, and provided that x is 
not too close to one of the classical turning points. It will also be convenient to 
write these functions in modulus-argument form: 

9± = V7^e'*± 

Finally we make the approximation 



1 



lS(x'±^]^ ^S{x) + y,,{x) (x'±^^x]+ i-p',,{x) [x'±'-~x 



1 



1 



y 



1 



2r 



This approximation will be justified provided 



T |p"ci(^)l X ) « 1 a^^d - \p\,{x)\ {mny « 1 



(26) 



(27) 



or, using the fact that |p"^j(x)| ~ m/(wa;,iiax) everywhere except in the vicinity of 
the classical turning points. 



mtob^ 



and Smuih b'^ <^ x„ 



(28) 



Making these approximations in Eq. ( P4[) and carrying out the Gaussian integrations 
gives, after a certain amount of algebra. 



W,:cd{t,X,p) « Wclit,X,p) + Wosc{t,X,p) 



(29) 



where 



Wci{t,x,p) ^ —^exp -a'i{p + pci{x)) p-{x) 



: exp 
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and 



Wosc{t,X,p) 



4fi.(Tl(T2V A 



exp 



-a2{p + PPciix)) - cri(pci(a;)) 
2 



X cos 



:S{x) +xit,x,p) 



p^{x)p+{x)Y (30) 



and where we have set 



P± = \9±\ 



4 



a ± 2cp^j(x) 
A 



h^aA + b{p'Jx)y 



h^aA + b{p'Jx)y 



h^a? + (1 - 2h?c'^ + h^A^) {p'^i{x)y 
c{l + h^A)p'^,{x) 
h?aA + b{p'^,{x)f 



h^bHp'jx)y 



x{t, x,p) is a phase whose functional form is unimportant for present purposes. 

Wci is non-negative, and it is concentrated on the classical energy surface at 
p = ±pci- It is therefore a possible classical phase space probability distribution 
describing a particle of energy En, with p+ (respectively p_) being the probability 
density function for the particle to be located at x and moving to the right (re- 
spectively left). On the other hand the term S/h in the argument of the cosine 
means that Wosc is very rapidly oscillating. Wosc is the term responsible for the 
tendency of the Wigner function to swing negative. It may therefore be regarded 
as the quantum mechanical correction to the classical distribution. 

It can be seen from Eq. (BO) that W^osc will become negligible once 



In that case W^ed 



criPciix) > 1 

Wci and, in view of Eq. (p3|), 

- , ^ Pci{x) ( p+{x) ~ p^{x) 
ve{x) « — 



TO \p+{x) + p-{x) ^ 

from which it follows that ue lies within the classical range 

Pc\{x) ^ _ , . ^ Pc\{x) 
< ve.{x) < 

TO TO 

Specialising to the case of the Caldeira-Leggett model it can be seen from 
Eqs. ([l9|), ( p5| ) and ( p8|) that the condition for the approximation of Eq. ( p6| ) to be 
valid is 



At 



(31) 



where is the time at which the smearing in the x-direction becomes significant 
[see Eq. (pT|)], tioc is the localisation time [see Eq. (p2[)], and Ab ~ h/{mLux^s^x), as 
before. We also have 



2 2 



1 



9to4?i2 
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Provided that x is not too close to one of the classical turning points dX x — ±a;n 
we have Pci{x) ^ mLox-a^a:^ and p' ^•^{x) ^ muj. Consequently 



2 2 



(if 



(a (A) 



1 + 4 



The fact that |'0svs) is highly excited means that a^max 3> Ab- Taking into account 
inequalities (|3|) we conclude that, in the case of the Caldeira-Leggett model, Wqsc 
is negligible, and we is approximately within the classical range of values, once 
t > ic- 

Finally, we remark that it follows from Eqs. and ( p^ ) that 

t\oc V ^^fi J 

where — {n + \/2)Tih} — {1 / 2)mLo'^ x1^^^ is the mean energy. We see from this 
that, in the case of a macroscopic body, is typically <^ tioc. As we mentioned 
above, tioc is the time at which the Wigner function becomes strictly non-negative, 
for every possible choice of initial state However, the argument just given 

shows that, for states of the form specified by Eq. (^), the Wigner function typically 
becomes approximately non-negative very much sooner than this, and approximate 
non-negativity is enough to ensure that we lies approximately within the classical 
range of values. 



4. Calculation of Wg"-* {t,x,xi, . . . ,xn) 

The requirement that IweI < Pc\/'m is a necessary condition for the Bohmian 
trajectories to be quasi-classical. However, it is clearly not sufficient. We therefore 
need to turn from the ensemble-averaged quantity v^{t, x) to the Bohmian velocity 
itself, Ug"-* {t,x,xi, . . . , xn). 

In order to calculate Wg"'' it is necessary to resolve the ambiguity mentioned in 
the Introduction, arising from the fact that the density matrix pbath [see Eq. (^], 
describing the initial state of the heat bath, does not uniquely determine a cor- 
responding ensemble [see the discussion in the paragraph following Eq. (^)]. One 
obvious choice is to represent pbath in terms of eigenstates of the heat bath Hamil- 
tonian: 

/3bath=5]pij|i?)(i?| (32) 
E 

where \E) is the eigenstate of ffbath with eigenvalue E and pe = 3\f exp [—E/{kBT)] 
[c.f. Eq. (|^)]. However, we will find it more convenient to use the coherent state 
representation 

Pbath = J dxidpi . ..dx^dpN P {xi,pi,. . .,xn,Pn) 

X . . .,xn,Pn) (Si, Pi, • • .^Xn^PnI (33) 

where P is the thermal Glauber-Sudarshan P-function p3^, R4l pa, pfl 



N 



[e^^ - 1) 



P{xi,pi,...,xn,Pn) = Y[ 1 exp 



r=l 



\2 
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and \xi,pi, . . . ,xn,Pn) is the coherent state with a;-space wave function 



(a;i, . . . ,XAr I a;i,_pi. 

N 



, xn,Pn) 



n 

r=l 



rrirUJr 



cxp 



rrirUJr 

2h 



2h 



[3r denotes the ratio hwr / {k-^T) ^ . 

Eqs. ( ^2| ) and (^3|) are completely equivalent ways of writing the density matrix. 
Consequently, from the point of view of the Copenhagen interpretation it makes 
no difference, whether we take the ensemble corresponding to pbath to consist of 
the states \E) with (discrete) probability distribution pe, or whether we take it to 
consist of the states \xi,pi,. . . ,X]S!,Pn) with (continuous) probability distribution 
P{xi,pi, . . . , xn tPn)- However, from the point of view of the Bohm interpretation 
it does make a difference; since in the former case, but not in the latter, the Bohmian 
velocity of the r"^ heat bath particle is certainly zero at t = 0. In the following 
we will assume that the ensemble giving rise to pbath consists of coherent states, 
since the integrals are then much easier to calculate. The question, as to whether 
choosing a different ensemble would require our conclusion to be modified, we defer 
to a future investigation. 

Eq. (|^), giving the state of system+heat bath at time t, thus becomes 



_ -itH/h 



sys/ 



Fi,Pi, 



where the label a now denotes the 27V-tuple (xi, 
sponding Wigner function: 



■ ,xn,Pn)) 

,Pn)- Let Wa be the corre- 



Wa(t,x,p,xi,pi, 
= ^ fdydyi.. 



■ ■ , xn,Pn) 
. dyN exp 



N 

E 

r=l 



Pry, 



2 ' ^1 



. . . , Xn 



^1 «'a(t))(*c.(i)| 2: + f,a;i + f , 



, Xn 



Vn \ 
2 / 

(34) 



and define 

Wa{t,X,p,Xi, 

Eqs. §), (H) and (|3 

„(") 



Xn) = j dpi. 
imply 

'{t,x,xi, . . .,xn) = 



■ dpN Wa{t,x,p,xi,pi, . . . ,xn,Pn) (35) 



/ dppWa{t,x,p,xi, . . .,Xn) 
m J dp Wa{t, x,p,xi, . . . , xn) 



(36) 



This equation is of a similar form to Eq. (p3|). It should, however, be noted that 
Wg"'' is the actual Bohmian velocity, whereas the quantity ue given by Eq. ( p3| ) is 
only an average. 

In order to calculate this quantity we follow Halliwell and Yu [Q, and use the 
fact that, because the Hamiltonian is quadratic in the positions and momenta, the 
Wigner function propagates in the same way as a classical phase space distribution. 

It is convenient to employ the vector notation: 



'n 



Vr 



Let ri(t), rij.{t) be solutions to the classical equations of motion, which result from 
the classical analogue of the Hamiltonian of Eq. (Ilj). Since the Hamiltonian is 
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quadratic we have 



N 



77(t) = A(<)r,(0) + ^B.(tK(0) 



r=l 
N 



r^M = a(t)r,(0) + D„,(t)77,,(0) 



(37) 



(38) 



for suitable matrices A(^), Br(t), Cr{i), T)rr'{i)- In Appendix ^ we give exact, 
closed form expressions for these matrices. In Appendix ^ we show how the matrix 
M appearing in Eq. ( p^ ) can be expressed in terms of them. 

We can use these matrices to propagate Wa forward in time |3^ : 

Wa{t,r},r]^, . . . ,r7jv) 

= ( 0, A(-t)r7 + Y Br(-i)»7r, Ci(-i)r/ + ^ Di,(-t)r/,, 



.,CN{-t)r] + Y^Nr{-t)r]A (39) 

r ' 



Also 



T/Fa(0,r7,?7i, . . . ,r7jv) = Wsys{'n)wlll^{'n^, ■■■,Vn) 



where Wgys, M^bath ^""^ ^^"^ Wigner functions corresponding to IV'sys), 1*71, 
respectively. We have [|2l|l|l 



■)N 



N 



where 



1 / rririOr 



Using Eqs. (||), (||), (|^), (|l]) we deduce 

Wait, r],xi, ... ,xn) ^ J (fr]'G{t,r],xi,...,XN\'n')Wsys{'n') 

where 

G{t, T],Xi, . . .,xnW) 

2N-2 



(40) 

(41) 
(42) 
(43) 



d $,dpi . . . dpN exp 



N 



N ^ 
- {^ri-t)v + T>rr'{-t)r]^, - fj^) Ar(Cri~t)r] + T>rr" i~t)r].^» ~ f), 

r,r' ,r" — 1 

(44) 

Carrying out the Gaussian integrations in Eq. (p4|), and using the results in Ap- 
pendix gives an exact, closed form expression for G. However, the expression is 
rather complicated, due to the fact that couples together the different oscilla- 
tors constituting the heat bath. We will therefore confine ourselves to the case when 
the interaction betweeen heat bath and system is weak. It is shown in Appendix |^ 
that we may then approximate 



A{t) 



cos ujt sin ujt 

'VrLLOsmLot cos Lot 



(45) 



13 



Brit) 



Crit) 



Kr ( mrhi°\t) hi°\t) \ 

mmrioujr ymiirirftr\t) mhi"\t)J 

Kr ( mhlP^t) h'l°\t) ' 

mrrirUJuJr \inmr'hi"\t) mrhi°\t) 



where 



and 



(46) 

(47) 
(48) 



Ur sin ut — u! sin u!rt 



Brit) 



COSWri 



■ sin ujrt 



-mrUJr sinUJrt COS LUrt 

If we also assume that <C 1, then we can further approximate 



A{t) « 1 



and 



UJt: 



Using these results in Eq. (EJ) and carrying out the integrations gives 



G(i, J7, xi, . . . , XN\'n') ~ const, x exp 



N 



r=l 



EnirljJr ( I \\'' 

—f—[Xr - qr[r])) 



X exp 



[rj' -ri-SfM~Ht) iv'-V-S) 



where 



(r,)= (D,(t)(f^,,-a(-t)r,))^ 



N 



(49) 

(50) 
(51) 



M{t) = 2h / duj' I{lo') 



i^uj' t— sin io't^^ 



f^Lj'i — sincj'ij ^1— cos uj't^ 



(^Lj't— sin uj't^ (^1 — cos uj't^ ^1 — cos u't^'^ 

'(52) 

and I{uj') is the spectral density function defined by Eq. (|^). Eqs. ( p3| ) and ( ^ ) 
then imply 



miLOi . . . m^LON 
T,N+2f^N detM(t)^ 



exp 



N 



—f—{^r-qr{v)) 



Wa {t,r],xi, . . . ,xn) 



where we have set <5 ~ 0, which will be justified if the coupling between system and 
heat bath is sufficiently weak. Define 



X /d'ry'exp -^{rj' - -n fM-^tXtj' ^ q) Wsysiq' 



a c 
c b 



and 



detM-i = A 



(53) 



We then have, by essentially the same argument as the one leading to Eq. (|29|), 
that if 



< Xn 



and 



9~ 3 

SniLuh ^ x„ 



(54) 
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then 

where 

W±{t,x,p,xi, . . . ,xn) 

= const. X exp 

and 

Wosc{t,X,p,Xi, . . .,Xn) 



Wc,^W-+W+ + Wo, 



(55) 



N 



const. X exp 



AT 



X^^T"^!^*- -^r{x,p)y ~ al{p + f3pci{x)y - di'^ {pci{x)y 



X COS 



-S'(x) +x{t,x,p) 



{P-ix)p+ix))- 



with 



p± = i5±r 



~a±2dp'^,ix) + b{p'^,{x)y 
A 

h^hA+b{p'^^{x)y 



h^hA + b{p'jx)y 



hH"^ + (i - 2ft2g2 + f^i^2\ (p'^^{x)y + h?b-^{p'^^{x)y 



c (l + h^K) p'^.ix) 



(56) 
(57) 

(58) 
(59) 

(60) 



n''aA + b{p'Jx)y 
x{t, x,p) is a phase. 

If (Tipci ^ 1, then T4^osc win be neghgible and 

For a given value of x the Gaussian in the expression for W± is peaked at the point 
p = ±pci{x), Xr — qr{x,±Pci{x)) . Supposc that Xr = qr{x,Pci{x)) for r = 1, . . . , N . 
Then 



W- — const. X exp 
where [c.f. Eq. @)] 



^2 



-o-l{p - Pc\{x)y - a^_{p + pci{x)y 



1 I Ti /-.{sin-Oj't — uj'tcosoj'tf 
aoj I[uj ) 



,/4 



(61) 



It foUows that, if (7_pci, o'^pd ^ 1, then will be negligible, so that 

Referring to Eq. ( |3^ ) we see that this imphes Wg"^ w +pci/m. Suppose, on the other 
hand, that Xr = qr{x, —pci{x)) for r = 1, . . . , iV. We then find that W+ will be 
negligible if a+pd, (r^pd 3> 1, in which case Wg"-* « —pd/m. The configuration space 
probability density function is obtained from Wa by integrating out the momentum. 
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It follows, that with probability close to 1, either Xr ~ qr{x,pci{x)) for r = 1, . . . , N , 
or Xr ~ q_r{x-, ~Pci(2;)) for r — 1, . . . , N. We conclude that there will be a high 
probability of being close to one of the two classical values ±Pci/m provided 

S-iPci > 1 cr±pci > 1 dspci > 1 (62) 

Let us now specialise to the case of the Caldeira-Leggett model, for which the 
spectral density defined by Eq. (|^) takes the form jl^ 

„ - ,- (63) 

10 il < LO 

for some cut-off frequency VL. Subsituting this expression in Eq. (^2|) gives |^ 



M(<) 




t{u~s\n li) (1— cos u) 
(1 — cos u) 



u 



7e) -;^(lnf^i-ln2 + 7E)\ , 



TT V (Inf7t-In2 + 7E) ' | Inf^t - i ln2 + |7e 

where 7e is Euler's constant. Referring to Eqs. ( p3[ ) we deduce that, if Infif ^ 1, 

a C\ TT 



c 



3m^ _ 

W * 1 and A 



2 



It can be seen that a, &, c and A are cut-off dependent (unlike the quantities a, b, c 
and A considered in the last Section). 

With these values inequalities (jsj) become 

{cjtf ( ' « iujf In nt « ( ^) ' (64) 

where Ab = h/ (mwxmax)- 

Eqs. iM) and (H) imply |37| 



X.2 



a- 



3 



47* r , (sinw - ucosm)^ 27i^lnf7t 

n;^h 



if In fit ^ 1 (where 7e is Euler's constant, as before). If x is not too close to one 
of the classical turning points Pc\{x) ^ mujXnia.x and p'cii^) ~ rnuj. Consequently 

~2 2 / \ 22:niax ,2 , oj- 

ttAr 



and [c.f. Eq. (|8|)] 

o-?P?i(a^) 



4a;niax CiJ7t^ In fit 



SttAb 1 + ^ (^7^2 In m)^ 

Taking into account the fact that these equations assume that t is in the range 
specified by inequalities (Q), we see that aipd and aapd will be ^ 1 provided 
Lojt"^ In fit > Ab / a;max • 

Referring to Eq. ( |57| ) we see that, away from the classical turning points, 

(since we are assuming tot ^ 1). Consequently a±pci ^ 1 if Inili ^ t^^a^max/ (t^b)- 
In the case of weak coupling (so that 7 <C w) and large quantum numbers (so that 
Ab ^ a^max) this inequality is automatically satisfied, for all physically reasonable 
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values of fit. The classicality conditions ( |62| ) then reduce to the single requirement 
t ^ ic, where tc is the solution to 

0-17^^ In nic — — — 

Finally, let us consider the relation between tc and the quantity tc discussed 
in the last section [see Eq. (pi])]. We clearly ought to have ic > tc] however, it 
is not immediately apparent that this is necessarily the case. In fact, the seeming 
difficulty disappears once it is recalled that, in deriving the Caldeira-Leggett master 
equation from a Brownian motion model, it is assumed [ p^ that fceT 3> hfl and 
V.t ^ 1. In particular, the discussion in the last Section tacitly assumes fltc ^ 1. 
Hence 

which implies that ic ^ ic- 



5. Conclusion 

The motive for this investigation was the question, whether it is true that any 
process which tends to produce decoherence also tends to make the Bohmian trajec- 
tory of a macroscopic object approximately classical. The results we have obtained 
provide some support for this hypothesis. However, it would clearly require more 
work to settle the question. We have only considered a particular, idealised model 
of the interaction between a macroscopic body and its environment. Moreover, 
our results were obtained on the assumption that the ensemble described by pbath 
consists of coherent states [see Eq.(|33|), and discussion in paragraphs following]. It 
would clearly be desirable to see if similar results hold in the case of other models, 
and for other choices of ensemble]^. 

However, realistic models of the interaction between a macroscopic object and its 
enviroment are very complicated, so that detailed calculations, of the kind carried 
out in this paper, are not usually feasible. What one needs is a general princi- 
ple, or mechanism, which can be shown to be operative even in those cases where 
the complexity of the problem makes detailed calculation impracticable. The most 
promising candidate for such a mechanism is the process of environmental moni- 
toring. 

Particularly relevant in this respect is a recent paper by Halliwell [^ . HalliwcU 
analyses Brownian motion models of the kind considered in this paper, and he shows 
that the positions and momenta of the heat bath oscillators constitute a store of 
information about the trajectory of the system particle. He also shows that there 
is a relationship between the amount of information stored in the environment and 
the amount of decoherence. One may plausibly speculate that a similar principle 
holds true with regard to the Bohm Interpretation: namely, that there is a direct 
relationship between the degree to which the Bohmian trajectory is approximately 
classical, and the amount of information about the trajectory which is stored in 
the environment. It would also be interesting to know whether such a principle 
applies to some of the other interpretations which have been proposed in which the 
particles follow determinate trajectories [|o[ |l[ |2[ |3|, Q. 



^ In this connection we should mention a recent paper by Geiger et al |38[ , in which the au- 
thors attempt to derive approximately classical Bohmian trajectories by making certain postulates 
regarding the form of the many-body wave function describing a macroscopic object. 
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Appendix A. Closed Form Expressions for A, B^, and D^r' 

The analysis in Section || is based on Halliwell and Yu However, Halliwell 
and Yu do not give explicit expressions for the matrices A, B^, and D^r'- The 
purpose of this appendix is to derive such expressions. Define 

^\ / ^— \ fO 



mojo (J / X—irirUJr y " 



Referring to Eqs. (37) and (38), and to the classical analogue of Eq. {m, we see 



that A, Br, Cr and Dj-r' satisfy 
d ^ d ^ 

— A{t) = LUQCrQA{t) — K.r<T_Cr(i) — Br(<) — UjQO-Q'Brit) — Kj.' (T -Gr' r(t) 

(Jib (JjL 

r—1 r—1 

subject to the initial conditions A(0) = 1, Br(0) — 0^(0) = and Drr'(O) = (5,.r'l- 
It is convenient to re-write these equations in integral form: 



(65) 
(66) 



A{t) = e'^°"'° -^nr f dt' e""'(*-*')'^«cr_Cr(i') 
^r{t) = -^Kr' / dt'e"''(*-*')'^°cr_Drv(t') 

r' = l 

Cr{t) = -Hrl dt'e'^'-(*-*')'^'-cr_A(t') (67) 

(68) 



T>rr'{t) = Srr'e'^^'"^ - Kr f dt' e'^^'^'-''^"^ a -Br' (f) 

Jo 

Eqs. (H) and imply 

A{t) = e'^°*'^° + / dt'Lit - t')A{t') 
Jo 

while Eqs. (|6^) and (|8|) give 

Br{t)^-Kr f (it'e"°(*-*')'^°cr_e'^'-*''^'- + /" dt'L(t - i')Br(i') 



(69) 



(70) 



where 

AT 



gC^0(t--t')<T0^_gC^.i'<T„^_ 



Br can be expressed in terms of A: 



Br{t) ^ ~Kr [ dt' A(t')(T^e' 

Jo 



LJr{t-t')iTr- 



(71) 



as can be verified by substituting this expression into Eq. ( |70D and using Eq. (|& 
Carrying out the integration in the expression for L wc find 



L(t) 



( x{t) 

\rnx{t) 
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X[t) = / duj I{uj ^2 



where 



mujo J„ 

being the spectral density defined by Eq. (||). Let g be the solution to the 
integral equation 

5(t) = sinwoi+ / dt' x{t-t')g{t') 



In terms of this function the solution to Eq. (k39) is 



A(t) = ^f"f?) (72) 



Eqs. (|67|), (|68|) and (|7lJ) then imply 

mrrirUjQLUr \mmrnr(t) mhryt) J 
mrrirUJoLUr \mmrnr(t) mrtiryt) J 

jy„'{t) = J^-e--*-- + ( ^r'frr'it) frr'{t) \ (75) 

mmrnir'OJoUJr^r' \mrmr' Jrr' [t) TnrJrr'{t)J 



where 



hrit) ^ - [ dt' g{t-t') sinujrt' 



f , , UJrSmLOr't' -UJr'SmUrt' 

Jrr'[t)= dt g{t-t) 2 2 

Jo - ^r' 

We see that all four matrices may be expressed in terms of the single function g. 

Eqs. ([72|-f75|) are exact. Let us now consider the case of weak coupling. Working 
to first order in the we have 

g{t) = sinwo^ + O(k^) 

and consequently 

A(<) = e"°*'^° + 0{k^) 

mnirUiQijJr XmrUrhr (t) mh). ' (t) I 

mmrLUQUJr ynrurhr {t) nirhr ' (t) I 



where 



h^r\t) = - / dt' sm.uJo{t-t')si-a.uJrt' = - 
Jo 



LOr sin UJot — LOo sin LOrt 



.2\ 



We also note that the frequency counterterm is 0(k ), so = ^ to this order of 
approximation. This proves Eqs. (p5|-|4^). 
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2. Expression for M in Terms of A and 

The purpose of this appendix is to derive the relationship between the matrices 
A, Br, Cr, Drr', and the matrix M which appears in the integrated form of the 
Master Equation, Eq. (16). We wiU also show that the matrix A appearing in 
Eq. (|l|) is the same as the matrix A derived in Appendix 

We begin by noting that D^r', regarded as a 2N x 2N matrix, is invertible. In 
fact, it follows from the time-reversibility of the classical equations of motion that 

N 



Ait)A{-t) + J2 Br{t)Cri-t) ^ 1 (76) 

r=l 

N 

Cr(t)B,./(-t) + 'Drr"it)'Dr"r'i-t) = (5„/ 1 (77) 

r" = l 

N 

A{t)Bri-t) + Br/(t)Drv(-i) ^ (78) 

r' = l 
N 

Crit)A{-t) + J2 Drr' (0^,-' (-t) = (79) 



r' = l 

It is then straightforward to verify that 

N 

Y,'Drr"{t)-D;,l,{t)^Srr'l (80) 
r" = l 

where 

D-l(i) = Brr'ht) - Cr{-t)A-\-t)Br'{-t) (81) 

As before, we assume that at time t system+environment are described by the 
density matrix 

where pbath is the thermal state defined by Eq. (||). Let W{t, ri,rii, . . . rj]^) be the 
corresponding Wigner function. Then 

Wrcd(i, 'n) = J ■ ■ ■ d^VN Wpit, J7, J7i, . . . , J7^) 

We have, by essentially the same argument as the one leading to Eq. (^3[), 

Wredit,v) = J dVG(t,r7|r?')Wred(0,V) 

where 
Git,v\v') 



N / „ N 



^■X I d?id^r^^ . . . (fr^j^ exp - ^ Uanh(^) (pr{-t)'n + ^ -Drr-{-t)-n. 

J V r=l ^ r' = l 

Af s N ■ 

Ar(Cr{-t)ri+ Y D„,(-t)r;,,) J +^C^(^7'- A(-Oi7-^B,(-Oi7r) 

(82) 



r=l 



and where A^ is the matrix defined by Eq. (p2|), f3r denotes the ratio hujr / {kBT) 
and 3C is a normalisation constant. 
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Making the substitution ry'^ = Cr{—t) + J2r'=i'^rr'i—t)ri^, in the integral on 
the right hand side of Eq. (plh we obtain 



X exp 



G{t, r]\T]') = const. / d^^d^r]\ . . . d^r/ ^ 

^(^tanh(^)r,'jA,V, - ^eA-\t)Br{tW 



N 

El 



where we have used Eqs. (8^) and (^T]), together with the relations [which are easily 
seen to follow from Eqs. (T^-^l])] 



N 

E 

r' = l 



N 



Ai-t)~ ^ B,„H)D;,i(-t)a(-t) = A-i(t) 

r.r' — l 

Carrying out the Gaussian integrations we deduce 



G{t,ri\ri') 



1 



det AvdetM 



exp 



77'-A-i(t)r,) M-\t)[r,' - A-\t)'n 



(83) 



where 

m=A-\t) (^f]coth(^)B.(i)A,Ti(B.(t))^^ {A-\t)f 
and where the normalisation constant is fixed by the requirement 

Comparing with Eq. ( p^ ) we see that the matrix M appearing in Eq. ( p^ is the 
same as the matrix given by Eq. (§1). 



[9: 
[lo: 
[11 

[12 
[13 

[14 

[15 
[16 
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